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Abstract 


This  progress  report  deals  with  various  aspects  of  the  experimental 
work  and  the  theoretical  work  performed  under  this  contract. 

In  the  introduction  we  summarize  the  objectives  of  the  contract. 

In  the  experimental  work  we  first  describe  the  two  noise-measuring 
systems  which  were  employed.  A  Hewlett-Packard  digital  fast  Fourier  trans¬ 
form  spectral  analyser  was  purchased  and  adapted  for  our  purpose  by  writing 
several  computer  programs  for  the  output  display.  Test  jigs  for  devices 
under  test  and  preamplifiers  were  built.  Further,  a  real  time  dual  channel 
correlation  spectral  analyser  was  designed  and  partially  built. 

Measurements  involving  temperature  fluctuations  and  correlations  in 
metal  films  and  in  integrated  circuit  transistors  are  described.  The  latter 
have  been  recently  completed.  It  is  conclusively  shown  that  temperature 
fluctuations  are  not  the  cause  of  1/f  noise  in  these  devices. 

Under  the  heading  of  theoretical  work  an  extensive  analysis  is  pre¬ 
sented  of  noise  in  diffusion  systems.  Several  analytical  models,  not  con¬ 
sidered  before,  have  been  solved  and  the  noise  is  presented.  It  is  con¬ 
clusively  shown  that  volume  (bulk)  diffusion  noise,  whether  based  on  the 
physical  diffusion  source  or  on  Voss  and  Clarke's  P-source,  never  leads  to 
1/f-type  spectra.  However,  surface  fluctuations  as,  e.g.,  caused  by 
stochastic  surface  generation-recombination  processes,  can  give  rise  to 
many  decades  of  1/f  noise.  This  is  worked  out  in  detail  for  a  MOSFET 
structure.  We  conclude,  therefore,  that  the  large  noise  observed  in  MOSFETs 
can  be  due  to  carrier  number  fluctuations  subject  to  diffusion;  at  low 
frequencies  the  MacWhorter-van  der  Ziel  model  may  be  operative,  whereas  at 
higher  frequencies  our  model  takes  over. 

Finally,  in  the  theory  section  we  present  a  novel  derivation  of 
quantum  1/f  noise,  occurring  in  emission  processes. 


1.  Transistor  test  circuit  for  low-frequency  digital  spectral  analyser. 

2.  One-third  octave  display  of  base  channel  noise  for  transistor  2N3945  // 1. 

3.  One-third  octave  display  of  collector  channel  noise  for  transistor 
2N3945  #1. 

4.  Block  diagram  for  real  time  correlation  measurement  system. 

5.  Two-stage  FET  low-noise  amplifier. 

5A.  Active  load  low-noise  preamplifier. 

6.  First  main  amplifier. 

6A.  Second  main  amplifier. 

7.  Layout  of  films  for  correlation  measurements. 

8.  Noise  spectrum  of  aluminum  film  at  1.171  mA. 

9.  Noise  spectrum  of  aluminum  film  at  1.071  mA. 

10.  Noise  spectrum  of  aluminum  film  at  260  yA. 

11.  Top  view  of  integrated  circuit  CA  3018. 

12  (a).  Circuit  for  thermal  transfer  function. 

(b) .  Circuit  to  obtain  noise  spectra  and  coherence. 

13.  Thermal  transfer  function. 

14.  Correlation  of  RCA  3018  NPN  transistors. 

15.  Three-transistor  heat-transfer  measurement  setup. 

16.  Thermal  transfer  function  for  the  two  "sensors." 

17.  Correlation  (coherence  factor)  of  the  two  "sensors." 

18.  Noise  of  the  two  "sensor"  transistors. 

19.  Correlation  (coherence  factor)  of  the  two  "sensor"  transistors. 

20.  Low-frequency  one-third  octave  display  of  base  channel  noise  for 
2N3945  ifl  (see  also  Fig.  2). 

21.  Low-frequency  one-third  octave  display  of  collector  channel  noise 
for  2N3945  #1  (see  also  Fig.  3). 

22.  Low-frequency  coherence  between  collector  channel  and  base  channel 
noise  for  2N3945  #1. 


23.  Intermediate  frequency  coherence  between  collector  channel  and  base 
channel  noise  for  2N3945  it  1. 

24.  Noise  for  embedded  sphere;  no  volume  sources  n  4  0. 

25.  Noise  for  nonembedded  sphere  n  =  0. 

26.  Noise  for  nonembedded  sphere  n  /  0. 

27.  Noise  spectra  for  nonsymmetrical  bar  based  on  Voss  and  Clarke's 
P-source . 

28.  Noise  spectra  for  nonsymmetrical  bar  based  on  the  physical  diffusion 
source . 

29.  Noise  spectra  for  disc  and  square  bar  based  on  the  physical  diffusion 
source . 

30.  Noise  spectra  for  a  MOSFET  due  to  stochastic  surface  generation- 
recombination  processes. 
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A.  Introduction 

Whereas  many  noise  phenomena,  like  shot  noise,  thermal  noise  and  genera¬ 
tion  recombination  noise,  are  well  understood,  1/f  noise  remains  an  enigma. 
This  noise  has  been  observed  in  all  semiconductors,  metal  films,  and  semi¬ 
conductor  devices.  Besides,  this  noise  has  been  found  in  electrolyses  and 
membranes,  in  biomedical  phenomena,  etc.  Though  some  investigators  believe 
that  there  must  be  some  universal  phenomenon  operative  in  all  these  manifes¬ 
tations,  like  scale  invariance  (S.  Macklup  [1]),  we  believe  that  the  ultimate 
explanation  will  require  the  involvement  of  specific  noise  causes,  as  for  any 
other  type  of  noise.  This  does  not  mean  that  general  principles,  like  scale 
invariance,  or  nonlinear  coupling  (Handel  et  al.  [2])  are  irrelevant.  How¬ 
ever,  there  must  be  physical  mechanisms  underlying  these  mathematical  descrip¬ 
tions. 

To  date  there  are  four  serious  contenders  in  the  theoretical  arena. 
First,  there  are  still  many  attempts  to  explain  1/f  noise  in  terms  of  diffu- 
sion-^/ke  transport  phenomena.  In  particular,  the  possibility  of  temperature 
fluctuations  as  the  cause  of  noise  in  metal  films  has  been  often  suggested 
in  the  recent  past,  notably  by  Voss  and  Clarke  in  1976  [3]  and  later  by  the 
Chicago  group  of  Hover  et  al.  [4-6].  The  possibility  of  obtaining  1/f 
noise  from  surface  initiated  heat  transfer  was  examined  prior  to  this  grant 
by  van  Vliet,  van  der  Ziel  and  Schmidt  [7].  A  very  complete  further  investi¬ 
gation  of  the  possible  connection  of  diffusion  noise  and  1/f  noise  was  made 
under  this  grant  by  Mehta  and  van  Vliet  [8];  see  Part  C  of  this  report. 

Second,  there  are  still  attempts  to  explain  1/f  noise  in  terms  of 
number  fluctuations,  e.g.,  due  to  tunneling  to  deep  traps  and  "slow  surface 
states,"  as  in  the  van  der  Ziel  -  McWhorter  theory  [9],  In  this  case  the 


scale  invariance  of  1/f  noise  stems  from  the  very  large  distribution  in 
lifetimes  with  probability  P(t)  -  constant  /l.  This  type  of  theory  might 


be  applicable  to  1/f  noise  in  MOSFETs,  together  with  the  "fast  surface 
states"  mechanism  suggested  by  Mehta  and  van  Vliet  [ 10 1 . 

Third,  there  are  many  indications  to  date  that  in  many  devices  (per¬ 
haps  MOSFETs  excepted)  the  1/f  noise  is  due  to  mobility  fluctuations.  This 
was  first  suggested  by  the  pioneering  work  of  Hooge  and  Gaal  [11],  Hooge 
and  van  Damme  [12],  and  Kleinpenning  [13].  So  far  only  one  reasonably 
satisfactory  theory  has  been  presented  to  explain  mobility  fluctuations  in 
terms  of  phonon  number  fluctuations,  see  Jindal  and  van  der  Ziel  [14]. 

Fourth,  there  is  a  fundamental  1/f  noise  theory,  which  attributes  the 
noise  to  quantum  (infrared  divergence)  effects  in  scattering,  involving 
either  bremsstrahlung  or  the  excitation  of  other  elementary  excitations  and 
configurational  states  or  correlated  states  (Handel,  Ngai  [15,  16]).  This 
theory  also  interprets  the  noise  as  mobility  fluctuations.  Many  details 
are,  however,  poorly  understood;  in  particular,  the  magnitude  of  the  noise 
is  not  easily  evaluated  quantitatively. 

Under  the  present  grant,  we  have  started  with  the  examination  of  the 
first  hypothesis  which  involves  temperature  noise.  In  part  B.II.  we 
describe  experiments  which  investigate  the  possibility  of  temperature 
fluctuations  as  being  the  cause  of  1/f  noise  in  semiconductor  devices  and 
films.  Though  a  full  investigation  in  the  entire  temperature  range  from 
4K  to  300K  is  still  underway,  the  experiments  to  date  indicate  clearly  that 
at  room  temperature  the  temperature  noise  hypothesis  fails.  Similar  con¬ 
clusions  have  at  the  same  time  been  reached  by  two  other  research  groups  at 
Cornell  (W.  Webb  et  al.  [17]  and  at  the  University  of  Illinois  (Weissmann 
[18]).  The  concurrently  made  theoretical  studies  are  reported  in  Part  C.II. 

The  quest  for  mobility  fluctuations  will  be  aided  by  the  study  of 
1/f  noise  in  near-ballistic  submicron  devices,  obtained  from  Dr.  L.  Eastman 
at  Cornell.  Under  the  present  grant  the  amplifiers  and  correlation  noise 
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equipment  for  this  purpose  were  constructed,  see  Part  B.I. 

The  quest  for  the  validity  of  Handel's  quantum  1/f  noise  theory  has 
begun,  and  the  theory  is  being  studied  presently  under  the  extension  of 
the  OSR  grant;  the  experiments  which  will  be  decisive  for  the  rise  or  fall 
of  this  type  of  theory  are  described  in  the  extension  proposal  submitted 
recently  to  OSR. 

We  note  that  though  the  grant  started  in  November,  1979,  actual 
measurements  were  not  commenced  until  spring,  1980,  due  to  the  fact  that 
the  noise  lab  was  shut  down  after  the  leaving  of  Dr.  Rucker  when  Professor 
Chenette  was  at  the  NSF  in  Washington.  As  a  consequence,  the  start  was 
slow,  as  reflected  also  by  the  low  rate  of  dispensing  funds.  Consequently, 
a  no-cost  extension  was  granted  for  seven  months,  bringing  the  contract 
period  to  a  close  on  May  31,  1981. 

B.  Experimental  Work 

I.  Construction  of  Spectral  Analyzers  (R.R.  Schmidt) 

For  transistor  noise  measurements  at  low  frequencies,  two  systems 
were  designed.  The  first  was  a  digital  system  using  a  Hewlett-Packard 
3582A  Spectrum  Analyzer  which  calculates  the  fast  Fourier  transform  of 
the  time  record  for  frequencies  from  20  mHz  to  25  kHz.  The  second  system 
was  designed  to  cover  the  frequency  range,  3  Hz  to  30  kHz,  using  the  real¬ 
time  correlation  method. 

a.  Digital  low-frequency  system 

A  transistor  bias  and  preamplifier  test  circuit  is  shown  in  Figure  1. 
The  base  of  the  DTJT  is  AC  open-circuited  and  the  collector  is  AC  short- 
circuited  by  the  use  of  operational  amplifiers  in  noninverting  and  inverting 
configurations,  respectively.  The  preamplifiers  are  Burr-Brown  0PA101CM 
operational  amplifiers  which  have  an  equivalent  noise  resistance  of  2K  ohms 
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above  the  1/ f  noise  region.  The  DUT  operating  point  can  be  controlled  by 
the  variable  potentiometer  which  determines  the  base  voltage.  To  achieve 
maximum  low-frequency  response,  coupling  capacitors  have  been  avoided. 

The  noise  spectra  of  a  2N3945  npn  BJT  are  shown  in  Figures  2  and  3. 

The  high  frequency  plots  include  a  flat  region  and  excess  noise  regions. 

These  measurements  are  further  described  in  detail  in  Section  B.III.  of  this 
report.  Suffice  it  here  to  say  that  later  on  better  representations  for 
displaying  the  1/f  noise  were  developed, 
b.  Real-time  correlation  system 

A  block  diagram  of  the  real-time  correlation  method  measurement 
system  is  shown  in  Figure  4.  There  are  two  parallel  amplification  channels, 
the  outputs  of  which  are  multiplied  using  a  four-quadrant  multiplier.  The 
result  is  averaged  and  read  out  on  a  DC  voltmeter.  This  system  operates  at 
a  single  frequency  determined  by  the  bandpass  filter.  It  is  desired  to 
build  a  parallel  array  of  these  systems  operating  at  various  frequencies 
from  a  few  Hertz  to  30  kiloHertz.  At  present  a  prototype  adjustable  single¬ 
frequency  system  has  been  constructed  using  Spencer-Kennedy  Laboratories, 

Inc.  model  302  dual-channel  variable  electronic  filters.  Since  these  filters 
accept  high-level  input  signals,  in  the  prototype  system  they  have  been 
moved  after  the  second  main  amplifier  and  before  the  multiplier  in  each 
channel.  This  also  reduces  the  effect  of  60  Hz  hum  at  the  output  of  the 
filters . 

The  two-stage  FET  preamp  is  shown  in  Figure  5.  The  FETs  are  n- 

channel  2N4393s  with  high  g  and  low  noise.  The  switch  S,  allows  the  noise 

m  1 

figure  of  the  preamp  to  be  accurately  measured  by  comparing  the  noise  at  the 
output  due  to  thermal  noise  at  the  input  of  various  resistors.  In  addition, 
there  is  a  direct  input  for  very  low-level  signals  of  magnitude  the  order 
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of  100  nanovolts.  The  power  supply  V  must  he  capable  of  delivering  20 
milliamperes  at  about  250  volts  with  very  stood  regulation.  Hum  and  noise 
introduced  at  this  state  is  greatly  amplified.  The  40-kohm  load  resistors 
dissipate  two  watts  each.  The  low-frequency  response  of  the  svstem  is 
limited  by  coupling  capacitors  in  this  and  other  stages. 

An  alternate  low-noise  preamplifier  is  show  in  Figure  5A.  The 
input  stage  is  a  2N5087  pnp  transistor  in  the  common  emitter  configuration 
with  an  active  load.  The  resulting  gain  is  about  70  db.  A  JFET  source 
follower  using  an  n  channel  2N4393  with  very  high  input  impedance  directly 
couples  the  amplified  signal  to  an  output  stage  consisting  of  an  emitter 
follower  which  uses  a  high  6  npn  2N5088.  Only  the  JFET  stage  is  resistively 
loaded.  This  preamp  has  much  lower  power  consumption  than  the  previous 
preamp . 

The  first  main  amplifier  is  shown  in  Figure  6.  It  uses  a  Burr  Brow 
model  OPA102CM  low-noise  FET  input  stage  operational  amplifier  in  non¬ 
inverting  configuration.  This  results  in  a  high-input  impedance  that  does 
not  load  the  preamplifier  output  and  thus  yields  higher  gain.  The  gain  is 
set  to  about  60  decibels.  The  bandpass  filter  in  the  final  system  will 
be  a  low-level  active  filter  using  operational  amplifiers  or  discrete 
devices.  The  variable  attenuators  in  the  prototype  system  are  Hewlett- 
Packard  model  350Ds  which  will  be  replaced  by  smaller  built-in  units  in 
the  final  system. 

The  second  main  amplifier  is  similar  to  the  first  except  for  a  few 
differences  as  shown  in  Figure  6A.  There  is  a  600-ohm  resistor  at  the 
input  to  match  the  output  impedance  of  the  attenuator.  A  741-type  opera¬ 
tional  amplifier  may  be  used  since  the  signal  levels  in  this  case  are  much 
greater  than  the  background  noise. 
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The  multiplier  and  detector  used  in  this  system  may  be  the  same  as  is 
used  in  the  high-frequency  correlation  system. 

Under  the  extended  contract  these  dual-channel  analyzers  are  being 
completed. 


II.  Temperature  Fluctuations  and  Correlation  in  Kj.].raer  &  J  J  Lin) 

Thin  Films  and  Integrated  Circuit  Transistors 

a.  Noise  in  aluminum  films 


Aluminum  films  of  the  configuration  of  Fig.  7  were  prepared.  The 

0 


thickness  of  the  film  was  approximately  20  A.  With  the  other  dimensions 
as  in  Fig.  7,  and  a  specific  resistivity  of  2.5  x  10  ^S7cm,  one  finds  that 
the  resistance  is  of  the  order  of  4  k  ft,  as  was  observed.  Some  typical  noise 
spectra  of  one  of  the  films  are  given  in  Figs.  8-10.  For  Fig.  8 
(I  -  1171 u  A)  we  find  at  16  hz  with  Af  =  11  hz : 


sv  = 


(10-100.3/20)2 
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=  8.4  x  10-12  V2/hz 


With  a  mobility  of  300  cm2/V  ,  we  find  that  the  density  n„  and  the  number 

sec  y  0 

21  -3  10 

of  carriers  in  the  sample  are:  n^  =  8  x  10  cm  ,  -  7.5  x  10 

Thus  Hooge's  constant  follows  from 


SV  = 


otVz 
f  N 


or 


a  = 


SvfN 


We  find  a  -  0.46.  This  is  two  orders  of  magnitude  higher  than  the  commonly 
accepted  value  of  a.  We  therefore  attribute  the  noise  to  the  silver  paste 
contacts  to  these  very  thin  films.  Alternately,  the  noise  may  be  due  to  an 
oxide  layer  which  often  is  present  on  aluminum  films. 

The  correlation  between  the  two  films  proved  to  be  negligible  over 
the  entire  frequency  range.  This  is  what  would  be  expected  if  the  noise  is 
due  to  the  contacts  or  due  to  an  oxide  layer  (McWhorter  model) . 
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Similar  noise  measurements  have  been  made  on  tantalum  nitride  films. 
The  spectra  were  1/f.  No  correlation  measurements  on  these  films  have  been 
made  as  yet  due  to  technical  problems  in  producing  closely  spaced  films. 

b.  Noise^  in _thermallv _coupled _inte^rated_ circuit 
transistqrsj  _two  _transistor_  experiments^ 

Voss  and  Clarke  [3]  conducted  correlation  experiments  between  two 
adjacent  parts  of  a  metal  film.  From  the  observed  large  correlation  they 
concluded  that  the  noise  was  in  agreement  with  the  temperature  noise 
hypothesis.  However,  if  the  noise  is  really  due  to  temperature  fluctua¬ 
tions,  one  should  be  able  to  observe  a  large  positive  correlation  between 
thermally  coupled  but  electrically  isolated  parts  of  a  test  object. 

To  test  the  temperature  hypothesis,  we  examined  the  noise  and  thermal 

response  of  closely  spaced  transistors  of  several  RCA  CA3018  general  purpose 

integrated  circuit  chips.  One  chip  contains,  among  others,  two  separate 

transistors  and  a  Darlington  pair.  The  centers  of  the  separate  transistor 

circuits  are  about  200  pm  removed  from  the  Darlington  circuit;  the  mutual 

distance  of  the  separate  transistors  is  about  20  pm,  see  Fig.  11.  First  we 

examined  the  mutual  thermal  response  of  the  two  separate  transistors  by 

using  one  transistor  as  a  "heater"  and  the  other  one  as  a  "sensor,"  see 

Fig.  12. (a).  A  sinusordally  varying  signal  was  fed  into  the  base  of  one 

transistor  and  of  the  other  transistor  was  monitored;  according  to  the 
BE 

manufacturer,  V  is  a  sensitive  function  of  temperature.  When  a  heater 
BE 

signal  of  frequency  fg  was  fed  into  one  transistor,  the  other  transistor 
showed  an  a.c.  voltage  v^  of  frequency  fg  and  2fg.  The  second  harmonic  is 
a  consequence  of  the  quadratic  heater  power  Pg  <  cos2  2irfgt>, fed  into  the 
heater.  A  typical  thermal  transfer  function  is  given  in  Fig.  13.  The 
turnover  is  not  sharp  since  the  thermal  response  time 

Tth 


=  L2/a  , 


(1) 
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where  a  =  a  ,  /pc  is  the  thermal  diffusion  constant  (o  ,  is  the  thermal  con- 
th  th 

ductivity,  p  the  density,  and  c  the  specific  heat),  is  distributed  due  to 
the  variation  in  L  for  different  portions  of  the  device.  Assuming  the 
most  characteristic  to  be  20  pm,  taking  p  =  2.3,  c  =  0.7  Joule  /g°c,  and 

_  O  _  O 

oth  =  8  x  10  watt  /  cm°c,  we  find  a  =  4.8  x  10  cm2/sec,  t  ^  =  0.83  x  10 
sec.  This  gives  a  corner  frequency  of 


fQ  =  1/2lrTth  =  192  hz  ,  (2) 

which  is  of  the  right  order  of  magnitude.  The  shape  of  the  curve  of  Fig.  13 
can,  in  principle,  be  found  by  solving  the  diffusion  equation;  the  details 
of  such  a  calculation  will  appear  in  the  M.Sc.  thesis  of  J.  Kilmer. 

Next  the  noise  spectrum  was  measured  for  the  "sensor"  transistor,  in 
quiescent  condition  (no  heater  voltage  to  other  transistors) ,  using  the 
setup  of  Fig.  12.  (b) .  The  spectrum  showed  the  following  ranges: 


1/f  for  f  <  60  hz 
1/f*'2  for  60  <  f  <  1000  hz 
1/f  for  1000  <  f  <  10,000  hz 
white  for  f  >  10,000  hz 


The  spectrum  of  the  other  transistor  (the  "heater,"  but  without  heater 
power)  was  nearly  the  same.  Why  the  spectra  are  not  1/f  over  the  full  fre- 
quency  range  from  1  hz  to  10,000  hz  is  not  understood.  The  1/f 2  range 
indicates  the  possible  presence  of  diffusion  noise. 

However,  the  noise  is  not  likely  due  to  temperature  diffusion.  An 
order  of  magnitude  calculation  of  the  expected  noise  j^f  the  noise  was 
caused  by  thermal  coupling  to  temperature  fluctuations  of  the  chip,  using 
the  calibrated  data  of  the  thermal  transfer  function  in  Fig.  13,  indicates 
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a  much  too  low  noise  compared  to  the  observed  noise.  The  details  of  this 
calculation  will  also  appear  in  the  M.Sc.  thesis  of  J.  Kilmer.  That  the 
noise  is  not  due  to  temperature  fluctuations  was  further  substantiated  by 
measuring  the  cross  correlation  between  the  noise  of  the  two  transistors, 
see  Fig.  14.  The  coherence  factor  p  is  defined  as 

p  =  <v.v2>  /  \J  <v12>  <v22>  (3) 

As  is  noted,  p<0.1  for  most  frequencies. 

c.  Noise  in  thermally  coupled  integrated  circuit 
_transistors;  three  transistor  experiments 

The  two  transistor  experiments  have  the  drawback  that  a  strict  con¬ 
clusion  about  the  absence  of  temperature  fluctuations  as  cause  for  the 
observed  noise  is  not  possible  unless  a  quantitative  estimate  is  made  re¬ 
garding  the  expected  contribution  to  the  noise  and  to  the  correct  time 
coefficient,  based  on  the  calibrated  thermal  transfer  function.  To  cir¬ 
cumvent  this  problem,  we  devised  a  three-transistor  experiment.  In  this 
experiment  the  Darlington  pair  was  used  as  a  "heater";  the  collector  voltages 
of  the  two  other  transistors,  acting  as  "sensors,"  were  monitored,  see 
Fig.  15.  The  transfer  function  for  each  sensor  with  respect  to  the  Darling¬ 
ton  pair  heater  is  given  in  Fig.  16.  The  turnover  frequency  of  the  response 
is  now  of  the  order  of  2hz,  in  agreement  with  the  fact  that  the  transfer 
distance  L  =  200  ym  for  this  experiment. 

Next  we  measured  the  correlation  between  the  responses  of  the  two 
sensors.  The  result  is  given  in  Fig.  17.  We  notice  that  the  coherence 
factor  p  equals  strictly  1  for  frequencies  below  2hz;  above  this  frequency 
the  phases  of  the  two  responses  differ  due  to  different  thermal  delays. 

Finally  we  measured  the  noise  and  its  correlation  for  the  two  sensors 
under  quiescent  conditions  (no  signal  to  the  Darlington  heater).  The  spectra 
are  of  the  1/f  form  below  100  hz.  see  Fig.  18.  (At  the  time  of  writing 


- -  -  i  tfn  iHjJBii'dMiH  till  I  I'll  ii'il 
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the  noise  measurements  have  been  extended  down  to  frequencies  of  50  mhz, 
the  spectra  still  being  1/f.)  Lastly,  the  coherence  factor  for  the  noise 
of  the  two  separate  transistors,  used  as  sensors  before,  is  given  in 
Fig.  19.  Clearly,  p<0.02. 

The  coherence  of  the  noise  as  given  in  Fig.  19  can  now  be  directly 
compared  with  the  coherence  obtained  by  externally  induced  temperature 
fluctuations  as  found  in  Fig.  14 .  Clearly,  the  coherence  must  be  the  same 
for  all  temperature  fluctuations,  whether  externally  induced  or  spon¬ 
taneously  occurring,  since  the  small  signal  parameters  for  the  response 
are  the  same.  Since  p(Fig.  19)  <<p(Fig.  14),  we  conclude  that  the  noise 
cannot  be  due  to  spontaneous  temperature  fluctuations  of  the  integrated 
circuit  or  of  its  thermally  connected  parts.  Thus  temperature  noise  can 
be  excluded  as  a  cause  for  the  1/f  noise. 

The  three-device  experiment  is  clearly  a  very  powerful  tool  to 
investigate  temperature  fluctuations.  The  experiments  are  presently  being 
repeated  for  three  thermally  coupled  gold  films  of  submicron  dimensions, 
manufactured  for  us  by  Dr.  E.  Wolf  of  the  Cornell  Submicron  facility.  These 
experiments  will  then  be  done  for  all  temperatures  from  4  K  -  300  K.  This 
should  settle  the  uncertainties  about  the  effect  of  temperature  fluctuations 
as  possible  cause  for  1/f  noise. 

III.  Base-Collector  Noise  Correlations  in  Transistors  (R.R.  Schmidt) 

To  further  test  the  usefulness  of  the  possibility  for  measuring  the 
coherence  factor  of  two  noise  inputs  with  the  Hewlett-Packard  digital 
analyzer,  we  measured  the  base  noise,  collector  noise,  and  their  coherence 
for  a  2N3945  common  emitter  bipolar  transistor.  The  results  of  the  single¬ 
channel  noise  spectra  are  shown  in  Figs.  20  and  21.  We  note  that  there  are 

-2/3 

two  low-frequency  noise  regions.  Above  1  hz  the  excess  noise  goes  as  f  , 
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both  for  the  base  and  collector  noise.  Below  1  hz  the  base  noise  rises 
-3/2 

faster,  viz.,  as  f  ,  whereas  the  collector  noise  levels  off. 

The  coherence  factor  for  these  two  noise  spectra  is  given  in  Figs. 

22  and  23.  We  note  that  the  coherence  below  1  hz  is  rapidly  going  to 
zero,  whereas  above  1  hz  the  correlation  is  nearly  unity,  diminishing  again 
at  high  frequencies  since  white  noise  then  comes  into  play. 

The  coherence  factor  indicates  that  entirely  different  noise  mechan¬ 
isms  operate  below  1  hz  and  above  1  hz.  Thus  the  measurement  of  the  co¬ 
herence  factor  is  a  powerful  tool  to  separate  noise  phenomena,  to  find  the 
limits  of  the  range  of  a  given  noise  phenomenon.  This  is  particularly 
important  for  the  ubiquitous  1/f  noise,  of  which  it  is  usually  assumed  that 
no  lower  or  upper  bound  occurs.  In  the  transistors  shown  the  correlated 
1/f-type  noise  occurs  only  above  frequencies  of  1  hz,  up  to  the  white  noise 
regime  (several  kilohertz) . 

C.  Theoretical  Work 

I.  Transport  Noise  (K.M.  van  Vliet  and  H.  Mehta) 

a.  Stirvey_of^  .transport  H°is£ 

A  survey  of  transport  noise  methods,  involving  Green's  function  pro¬ 
cedures  and  eigenfunction  expansions,  is  appearing  in  Phvsica  Status 
Solidi  (b)  [19]. 

b.  Noi.s£  of_S}nnmetric£i._emb£ddec[  bod_i£s 

The  noise  of  symmetrical  embedded  bodies  can  be  found  in  closed  form 
using  the  Green's  function  procedure.  According  to  a  formula  of  van  Vliet 
and  Fassett  [20],  the  noise  spectrum  of  a  quantity  X,  subject  to  transport 
in  a  domain  Dg  of  volume  V  ,  is  given  by 

Sx(u>)  =  —  -yar  X  Re  /  /  G(r,iw,r')d3rd3r' 

s  D  D 

s  s 


(4) 
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where  G  is  the  Fourier-Laplace  transformed  Green's  function  of  the  trans¬ 
port  equation,  typically  the  diffusion  equation.  For  a  circular  or 
cylindrical  domain  the  solution  was  obtained  by  MacFarlane  [21],  Burgess 
[22],  and  bv  van  Vliet  and  Chenette  [23].  For  a  sphere  involving  diffu¬ 
sion  as  well  as  volume  generation-recombination,  the  spectra  were  hitherto 
unknown.  The  problem  was  solved  by  Mehta  and  van  Vliet  under  the  present 
contract . 

Thus,  we  consider  a  sphere  embedded  in  an  infinite  domain.  The 
diffusion  equation  reads 

4^-  +  ax  -  DV2x  =  0  ,  (5) 

dt 

where  a  =  t  ^  is  the  reciprocal  volume  lifetime.  The  Green's  function 
satisfies  the  Fourier-Laplace  transformed  equation 

-  ^  G(r,s,r ')  -  ^  G(r ,  s ,  r '  )  +  V2G(r,s,r')  =  -  ~  6(r-r')  (6) 


or 

(k2  +  V2)  G(r,s,r')  =  -  ~  <5(r-r' ) 

where 


For  a  spherical  coordinate  system  we  have 


72  . 

2  9r 

•• 


(r2  b  + 


Jf  r2sin2@  9cj>2  r2sin 


9  ,  .  r,  3  \ 

96  (sin  6  9?}  * 


(7) 

(8) 

(9) 


With  the  conditions  that  G  goes  as  for  R-*-o  (where  R  =  |r-r'|)  and 

K 


Urn  Gk(R)  =  0  , 

_  00  +  ico 
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one  obtains 


Gk(R)  =  e  1KR  ^4ttRD 


(10) 


To  find  the  spectrum  one  must  now  perform  the  integration  as  required  bv 
(4): 


1=1  f  d3rd3r'G  (R), 

V  V  K 

s  s 


where 

R  =  V  r2  +  r'2  -  2rr '  cos(8  -  8’) 


(11) 


(12) 


and 

d3r  =  r2  sin  GdQdcJ) 

d3r '  =  r'2  sin  8  ’  d0  *  d<t>  ’  . 


(13) 


One  may  choose  8*  =0.  After  much  algebra  the  integral  is  found  to  be 


4 

<Di 


2  2i<a  2i<a 

a^e  ae 


L  -2k^ 


i<; 


2i<a 


2k 


3iK 


2k2 


2k 


(14) 


where  a  is  the  radius  of  the  sphere.  For  the  spectrum  of  the  total  number 
of  holes  one  obtains,  with  X  =  P  =  /  xd3r: 

64 mp  a5  r  3 

S^p  =  Re  -  e  (<p2  -  4<J>  +  4)  +  (  +  <p2  -  4)  J  ,  (15) 


where  < j>  =  2i<a,  <  being  given  by  (8). 

We  first  considered  a  =  0  (pure  diffusion).  lie  result  is  then 
>5t 


AP 


16Trabp.  / 

=  -  \-2+y2+eY  [y2(sin  y + cos  y) 

Dy5  1 

+  4y  cos  y  +  2(cos  y  -  sin  y)  }  j  , 


(16) 
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where  y  =  2a  \  w/2D.  The  spectrum  is  plotted  in  Fig.  24,  see  the  n  =  0 
curve.  The  results  show  a  plateau  and  an  to  asymptote,  as  expected  from 
the  work  of  Lox  and  Orengert  [24].  For  diffusion  with  sinks  the  result 
is 


S4P<“> 


16ira5p0(g2  -  62)  +  16Tra5pQK 

3D(32  +  62)2  D(62  +  62)5 


(17) 


where 


K  =  (10B263  -  5g46  -  6s)  [266  +  26e  p  cos  6(6  +  2) 


-6 


+  e“6  sin  6(62  -  62  -  4  -  46)]  -  (65  -  106362  +  5B64) 

x  [62  -  62  +  26e~6  sin  6(6  +  2)  +  e"6  cos  6 (B2  -  62  +  4  -  46)  -4],  (18) 


and 


6  = 


6  = 


2  +  U)2 


2a 

ft 


=  [{^bT‘ 

=[  V-ia+i'(''2+“2  ] 


(19) 


(20) 


Furthermore,  the  following  dimensionless  quantities  were  used: 


u  =  lot,  ,  (t,  =  a2/D  is  the  diffusion  time), 
a  a 

n  =  t,/t  =  ax,  , 
d  v  d 


=  VT  V  h  +  \Tn* +  u2  , 


=  ^ ~2  ^  -n  +  ^  n2  +  u2  . 


(21a) 

(21b) 

(21c) 

(  2  Id ) 


The  spectra  appear  in  Fig.  22,  curves  for  n  ^  0. 

-2 

One  notices  the  occurrence  of  an  w  range,  in  addition  to  the  plateau 
-3 /2 

and  a)  range.  It  is  concluded  that  the  spectra  do  never  resemble  a  1/f 


spectrum. 
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c  •  Non  embedded^  sphere 

Next  we  considered  the  noise  in  a  sphere,  on  the  boundary  of  which  we 
have  surface  recombination.  The  boundary  condition  is  then 


D  —  G(r  =  a,s,r')  +  aG(r  =  a,s,r')  =  0 
an 


where  n  denotes  the  outer  normal  and  a  is  the  surface  recombination  velocity. 
The  differential  equation  for  the  Green's  function  is 


_2_  _3G  +  2  + 

r  3r 


1  r  i  3  /  •  s  % 

2  sin  6  39  (sin  6  36  ^ 
rz  L 


_i _  ' 

Ln26  3  <j>2 


where  k  is  the  same  as  in  Eq.  (8).  Using  an  expansion  in  spherical  poly¬ 
nomials,  the  solution  of  (23)  subject  to  (22)  can  be  shown  to  be 


where 


Qjl(<a)  = 


4irD  : 

u+1 

m| ) ! 

% 

,m 

[V" 

<> 

(xr') 

-  h 

(<a) 

-  (a  + 

DC. 

a 

DK^+1 

(tea) 

-  (o  + 

M 

a 

and  where  j  is  the  spherical  Bessel  function  of  order  l,  h  is  the 
spherical  Hankel  function  of  the  first  kind  of  order  l. 

We  must  now  perform  a  six-fold  integration,  as  required  by  (A).  The 
integral  is  obtainable  in  closed  form,  by  employing  the  addition  theorem 


-  21 


.  itcR 
l _ 

IkR 


00 


23  (2t  +  1)  j  (>cr) 

£  =  0  *  36 


(fr)  P  (cos  0) 

i  * 


(26) 


After  a  very  laborious  calculation  (see  the  Ph.D.  dissertation  of  Mehta  for 
details)  one  obtains 


Sp(u) 


4(var  P) 
D 


Re  'f'(cj)  , 


(27) 


where 


»(u)  = 


r)  +  iu 


1  - 


3j1(i V  n  +  in) 


C(n  +  iu)  j^(i  +  iu)  +  i-^n  +  iuj^i^n+iu) 


(28) 


where  n  and  u  are  defined  as  in  (21a)  and  (21b)  ,  while 


C  =  x  /  x  ,  t  =  a2/D  ,  t  -ala. 
s  d  d  s 

Here  t  is  the  surface  recombination  time  and  t,  is  the  diffusion  time, 
s  ci 

When  there  are  no  volume  sinks,  n  =>  0.  Using  the  relations  for 
spherical  Bessel  functions, 


jQ(z)  =  sin  z  /  z  , 

j^(z)  =  sin  z  /  z2  -  cos  z  /  z  , 

we  obtained 


Re  Y  (u) 


3  (AD  +  BC) 
4a4£(C2  +D2) 


where 


a  =  V  u  / 


2  and  where 


A  =  (-sin  a  cosh  a  +  a  cos  a  cosh  a  +  a  sin  a  sinh  a) 
B  =  (cos  a  sinh  a  +  a  sin  a  sinh  a  -a  cos  a  cosh  a) 


(29) 


(30) 


(31a) 

(31b) 
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C  =  (cos  a  sinh  a  +  a  sin  a  sinh  a -a  cos  a  cosh  a  -  co-— ^  s^-n--— —  )  (31c) 


„  ,  ,  ...  .  ,  sin  a  cosh  a  s  ... 

D  =  (-sin  a  cosh  a  +  a  cos  a  cosh  a  +  a  sm  a  sinh  a  H - - -  )  (31d) 


For  various  values  of  the  parameter  £  the  spectra  are  plotted  in  Fig.  25. 

In  contrast  to  the  embedded  sphere,  the  high  frequency  slope  is  always 

-2 

OJ 

When  both  volume  and  surface  g-r  processes  play  a  role,  q  ^  0,  £  /  0. 
The  following  result  is  found  after  considerable  algebra: 


Re  H’(w) 


£(MP  +  NQ)  -  3 (AM  +  BN) _ 

£  j  (?2P  +  12P)2  +  (B2Q  +  62Q)2  ( 


(32) 


where 


6  =  (1  /\TT)V  n  +  V  n2  +u2 

_  1 

(33) 

6  =  (1  /V  2)-y  -n  +  v  n2  +  u2  , 

(34) 

A  =  (-sin  6  cosh  B  +  6  cos  <5  cosh  B  +  3  sin  6  sinh  B)  , 

(35a) 

B  =  (cos  6  sinh  S  -  6  cos  6  cosh  8  +  6  sin  S  sinh  B)  , 

(35b) 

„  .  .  sin  6  cosh  8 

C  *  A  +  , 

(35c) 

_  _  cos  <5  sinh  8 

D  =  B  -  , 

(35d) 

P  *  (B2  -J2)  C  -  2bTd  , 

(35e) 

Q  =  (B2  -  B2)  D  +  2&I C  , 

(35f ) 

M  =  (B2  -S2)  P  -  2bTQ  , 

(35g) 

N  =  (B2  -  S2)  Q  +  2b  ?P  . 

(35h) 
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The  low-frequency  plateau  is 

lim  Re  f  (w)  =  (1  +  5£)  /  15 
ar*-o 

The  high-frequency  asymptote  is,  noticing  6  ~  6  '■\u/2, 
Re  T  (u)  3  /  C  u2  . 


(36) 


(37) 


The  possible  spectra  are  plotted  in  Fig.  26. 

Though  the  details  of  the  spectra  are  very  revealing,  we  note  that 
the  spectra  never  resemble  a  1/f  spectrum. 


d.  Non  symmetric  a_l_embedd_ed  bodies 

The  noise  was  also  calculated  for  a  nonsymmetrical  embedded  rectangular 
bar  of  dimensions  x  x  L^.  In  this  case  the  Green's  function  cannot 
be  obtained  in  closed  form,  and  the  eigenfunction  expansion  technique  is 
preferable.  Let  S  be  the  spectrum  of  the  volume  noise  source.  We  then 
obtain 


S(u) 
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(2tt)3 


/  /  d3k  d3k' 


,  ik. 
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(Dk2  +  io>)  (Dk'2  -  ico) 
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(38) 


For  Voss  and  Clarke's  P  source. 


sc(rrr2)  =  Pj,2^!-^) 


(39) 


For  the  rectangular  bar  we  then  obtain 


S(oi)  =  P^O^L^)2  / 


d3k  3  sin2(!Sk1L1) 


D2k4  +w2  i=l  (*5  k^  L^) 2 


oo 


(40) 
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For  the  bar  there  are  break  points  at  id.  =  D  /  L^2.  Voss  and  Clarke  claim 
that  the  following  frequency  ranges  occur 


w  <<  : 
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ot  ca 
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(40c ) 

a)  >> 
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cc  Cl) 

(40d) 

In  order  to  verify  this  Eq  (40)  was  rendered  in  a  dimensionless  form  and  put 
in  the  computer.  The  resulting  spectrum  is  given  in  Fig.  27.  We  note  that 
only  two  break  points  are  observable  and  that  n£  1/f  range  occurs. 

Voss  and  Clarke's  P  source  is  highly  hypothetical.  The  correct  source 
is  the  diffusion  source: 


S^(r,r’)  =  4D  V  •  V'  n(r)  6(r-  r') 


This  leads  to  the  result. 


S(w)  =  constant  / 

OO 


d  3kPk2 
D2kL|  +  w2 


3  sin2(%k.L.) 

n  - l-k- 

i=l  Ozk.L.)2 
1  1 


(41) 


(42) 


which  differs  from  (40)  by  the  factor  Dk2  in  the  numerator.  According  to 
this  formula,  the  following  ranges  occur: 


<<  0)^1 

:  S  (u) 

-  constant 

(43a) 

<<  O) 

<<  cu2 : 

S(io)  «  -  log  uj 

(43b) 

<<  0) 

<<  m3: 

,  ,  -1/2 

S(w)  «  u 

(43c) 

-3/2 

w  >>  S(w)  =  w 


(43d) 
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Again  the  spectrum  was  put  in  a  dimensionless  form  and  fed  into  the  com¬ 
puter.  The  result  is  given  in  Fig.  28.  The  predictions  of  Eqs  (43)  are 
more  or  less  confirmed. 

We  also  evaluated  the  spectra  for  a  square  bar,  <<  L^,  and  a 

square  disc,  >>  L^.  The  results  are  plotted  in  Fig.  29. 

Again,  we  conclude  from  all  these  computations  that  the  spectra  never 
show  a  1/f  resemblance  over  many  decades. 

In  conclusion:  We  established  beyond  doubt  that  diffusion  noise  gives 
no  1/f-like  noise  spectrum  in  any  geometry,  whether  embedded  or  nonembedded 
for  any  type  of  boundary  conditions  (Dirichlet ,  Neumann,  or  mixed).  Thus 
any  attempt  to  explain  1/f  noise  by  some  diffusion  mechanism,  with  or  without 
volume  g-r  processes  and  with  or  without  surface  recombination  is  doomed  to 
failure .  See,  however,  the  next  subsection  for  a  possible  mechanism  in 
MOSFETs . 

e.  In  the  above  calculations  we  considered  various  types  of  boundary 
conditions,  but  the  source  of  the  noise  was  due  to  diffusion  fluctuations 
and/or  volume  g-r  processes.  Next  we  consider  the  noise  due  to  surface  g-r 
processes.  In  this  case  we  have  stochastic  boundary  conditions,  i.e.,  (22) 
is  replaced  by. 


-  aAp 


+  aAp 


=n  = 


'  x=0 


'  x=w 


=  C(t) 


(44) 


Note  that  we  consider  one-dimensional  diffusion  in  a  sample  of  width  w.  The 
volume  source  £  will  be  set  =  0,  since  its  noise  was  calculated  in  the  pre¬ 
ceding  sections  (we  may  consider  the  noise  from  the  sources  £  -  surface  -  and 
£  -  volume  -  separately  since  these  sources  are  uncorrelated) .  The  one- 
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dimensional  diffusion  equation  for  holes  is 

3Ap(x,t)  _  D  32Ap(x,t)  +  Ap (x , t)  Q  ^ 

3T  3x2  T 

Associated  with  these  equations  we  define  the  Fourier-Laplace  transformed 
Green's  function  G(x,x',iio)  and  the  surface  Green's  function  H(x,0,i<o). 
For  the  noise  due  to  the  source  z,  we  obtain 


SAp(x,x',w)  =  H(x,0,iw)  H(x' ,0,  -  iw) 


(46) 


The  source  spectrum  is  white: 

=  4Maap0  /  A  , 


(47 


where  pQ  is  the  average  density  of  holes,  A  is  the  area,  o  is  the  surface 
recombination  velocity,  a  =  var  P  /  Pq  and  M  is  a  modulation  factor  >_  1. 

In  a  MOSFET  we  must  find  the  averaged  fluctuations  in  the  channel  of 
width  <<  W.  Thus, 


~~  f  H(x,0,iw)  dx  2 
1  0 


(48) 


For  the  surface  Green's  function  we  obtained 


H(x, O.im)  =  (sinh  ~  ^  cosh  ^)  , 


where 


y /  iw + 1  /  t 


and 


„  ,,  ,  Dz  .  .  ,  W  ^  2D  ,  W 

D  =  (1  +  -  )  sinh  —  +  —  cosh  —  . 

r  <j2L2  E  aL  L 


(49) 


(50) 


(51) 


We  write  W  =  W^  +  W2  where  W^  is  the  width  of  the  channel  and  W2  is  the 
width  of  the  bulk  material  of  the  MOSFET .  Carrying  out  the  integration 


indicated  in  (48)  we  find  S^—  (u))  01  |f(iu>)|  ,  where 


_L_ 

W, 


f(ico)  = 


W, 


W, 


,  '1  ,  D’L  ,  W2  COSh  ( L  )  ~  1 

tanh  I  I7  tanh  IT - w- - 

cosh  (— ) 


.  ,  V  D  .  n’L  ,  W2  .  D  ,  W1 

(tanh  —  +  — )  r— 7-  tanh  —y  +  1  +  —  tanh  — 
L  OL  DL '  L ’  OL  L 


(52) 


Here  primed  parameters  refer  to  the  bulk  and  unprimed  parameters  refer  to 
the  channel. 

At  very  low  frequencies  | W^/L |  <<  1  and  Iw^L  j  <<  1.  Expanding  tanh 
(W^/L)  =  W^/L,  etc.,  we  find  for  low  frequencies  a  constant  plateau. 

At  intermediate  frequencies  ]  /L  |  <<  1  but  |l.'2  /L  |  >>  1  since  W?  »  W^ . 
Expanding  the  appropriate  hyperbolic  functions,  we  now  find  for  this  range, 

j  f  (icj)|2  =  xd'  /  ts'2u  (53) 

where  x  , '  =  W2  /  d'  -  W  2  /  D  and  x  '  =  W  /  a  -  W „  /  o  . 
d  l  s  2 

Thus  a  large  1/f  range  occurs. 

For  high  frequencies  |w^/l|  >>  1  and  | /L ’ |  >>1.  One  easily  con- 
-2 

firms  that  there  is  an  w  asymptote.  This  is  confirmed  by  computer  calcula¬ 
tions,  see  Fig.  30.  For  realistic  values  of  W^  and  W^,  one  finds  that  the 
1/f  range,  which  starts  for  u)>l/x,  can  continue  up  to  many  gigaherz.  Thus, 
this  noise  which  is  due  to  fast  surface  states  continues  the  noise  spectrum 
due  to  the  slow  surface  states,  occurring  for  w<l/x.  This  is  the  cause 
that  1/f  noise  in  MOSFETs  can  occur  even  for  very  high  frequencies.  We  note 
that  the  results  obtained  stem  from  the  response  of  the  diffusion  equation 


to  surface  sources;  no  distribution  of  lifetimes  needs  to  be  invoked. 
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Having  investigated  many  diffusion  mechanisms,  we  conclude  that  1/f 
noise  seldom  occurs,  except  in  the  case  of  a  small  conducting  channel  on  a 
large  substrate;  the  1/f  range  stems  from  the  diffusion  and  subsequent 
recombination  in  the  bulk. 

The  details  of  this  theory  are  being  published. 


II .  On  the  Theory  of  Quantum  1/f  Noise 

(P.H.  Handel,  K.M.  van  Vliet,  and  A.  van  der  Ziel) 

Under  the  contract  the  above-named  persons  wrote  an  article  on 
"Superstatistical  Emission  Noise,"  now  accepted  for  publication  in  Physica  A. 
In  this  article  we  treated  various  emission  noise  problems  with  compound 
Poisson  statistics,  viz.:  noise  due  to  emission  centers,  cathodolumin- 
escence  noise,  photon  noise,  and  quantum  1/f  noise.  In  essence,  quantum 
1/f  noise  is  treated  here  as  a  wave  interaction  noise;  the  wave  packet  here 
being  of  a  corpuscular  rather  than  ondular  nature,  as  with  photon  noise. 

Thus,  consider  a  beam  of  emitted  particles  (electrons)  undergoing 
scattering  in  the  emission  process.  Let  m(t)  be  the  number  of  emitted  par¬ 
ticles  per  second.  The  distribution  for  the  emission  of  M  particles  in 
an  interval  (t,  t  +  0)  is  then  written  as 


W(M,0) 


{m( t) } 


m(t ' )dt ' 


M 


M. 


t+e 

-  f  m(t ' )dt ' 
e  v[m(t) ] 


(54) 


This  compound  distribution  takes  into  account  that  the  average  number 
emitted,  /£  m(t')dt'  is  subject  to  fluctuations  due  to  quantum  mechani¬ 
cal  uncertainty;  P[m(t)]  is  the  probability  density  functional;  may  also 

(m(t) } 

be  replaced  by  the  functional  integral  / Dm(t)  if  m  is  considered  to  be  a  con¬ 
tinuous  variable.  We  use  the  bar  for  averages  over  M  and  <>  for  averages 
involving  m.  For  the  variance  AM.2  we  find  from  (54)  in  a  straightforward 


fashion: 
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"V  -  VMe-D  -  He(Me  " 

_  t+e  t+e  t+e  t+e 

=  M  +  /  /  dt1  dt-  <m(t.)m(t-)>  -  /  m(t  )dt  /  m(t?)dt 

0  t  t  1  Z  t1  t  “ 

e  6 

=  0<m>  +  f  f  <Am(t  )Am(t-)>  dt  dt  (55) 

0  0  i 

The  first  term  is  the  shot  noise  term;  the  second  term  will  be  called 
"wave-interaction  noise"  or  'Intensity-fluctuation  noise"  (see  van  Vliet 
125]. 

By  MacDonald's  theorem  the  spectrum  S^(to)  is  related  to  AM^2,  according 

to  [noticing  that  0_1M„  =  ~  ft+Q  m(t')dt'  5  mQ(t)  is  the  "short-time 

G  0  t  0 

average"  of  m(t)]: 

°°  - 

S  (u>)  =  2u  /  d6  sin  co0  —  [02  Am  2  ] 

Am  q  do  0 

CO  _ 

=  2u  /  d0  sin  w6  [AM  2]  .  (56) 

0  39  9 

Using  (55)  this  yields  two  parts.  The  first  part  is 


S.1  (u>)  =  2w  /  d0  sin  w9<m>  =  2<m>  (57) 

Am  Q 

where  we  used  Cauchy's  principal  value  for  the  improper  integral.  Thus 
(57)  indicates  shot  noise. 

The  second  part  we  write  as 

°°  .  0  0 

S11  (w)  =  2w  <m>2  /  d0  sin  w9  v-r-  /  /  <Am(t  )Am(t,)  dt,dt»  (58) 

m  0  39  0  0  1  2 

<Am(t^)Am(t2)>  is  the  normalized  autocorrelation  function  of  the 


intensity-fluctuation  noise.  With  some  transformations  (see  the  appendix  of 
the  paper  for  details)  one  shows  that 


0  9*  *  0  u 

/  /  Am(t  )  Am(t„)dt  dt9  =  2  /  du  /  dv  <Am(v) Am(0) >  .  (59) 

00  1  1  “  00 


This  gives 

II  °°  6 

S  (w)  =  4u>  <m>2  /  d0  sin  u>@  I  dv  <Am(v)Am(0)>  .  (60) 

m  0  0 

Integrating  once  by  parts  and  using  the  Wiener-Khintchine  theorem,  we 

obtain 

CO 

S*1  (u)  =  4  <m>2  /  d0  cos  u0  <Am(0)Am(O)>  =  <m>2  S  (w)  ,  (61) 


where  S,  (u)  is  the  normalized  intensity-fluctuation  noise. 

Am 

*f* 

Let  now  T(x,t)  and  T  (x,t)  describe  the  local  particle  field  in  the 
usual  way  of  quantum  field  theory,  being  each  the  sum  of  a  non-scattered 

*t* 

beam  and  a  scattered  beam.  In  addition  we  assume  that  tp  is  an  analytic 


signal,  obtained  by  adding  to  half  of  the  signal  one-half  of  its  Hilbert 


transform.  Then  the  spectrum  S  j.  (f)  is  zero  for  negative  frequencies.  Thus 

4. 

the  local  density  (Jj'iJj  can  be  written  as  an  integral  over  positive  energy 

I 

i 

|  losses,  e  =  hf, 


<4/+vi/>  =  /  s  u(c)d£  ;  (62) 

0  >>' 


"f*  *!* 

hence  S  +(e)  is  both  the  spectral  density  of  tj;  and  the  part  of  <i^> ' H/>  with 
r 

average  energy  loss  in  (e,e+de),  divided  by  de.  For  the  latter  we  have, 


see  Handel  [15,16], 
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6(e)  +aA(e/en)aAe_1u(e  -  en) 

s  t(e)  - - - -2 - ^ - - 

1| <'  1  +  aA  /A  (e/e0)  e  xde 

e0 


where  aA  is  the  infrared  exponent,  u  is  the  unit  step  function,  and  =  hf^ 
is  an  energy  resolution  limit  set  (arbitrarily)  by  the  lowest  measured 
frequency  f  ^  T  ^ ,  T  being  the  duration  of  the  experiment;  A  is  the 
largest  possible  energy  loss,  i.e.,  the  kinetic  energy  mv2/2  of  the  carriers. 
The  delta  term  represents  the  "DC"  component  of  the  spectrum,  which  is 
therefore  affected  by  the  choice  of  e^.  Note  that  (62)  is  normalized  to 
unity: 


<i|ji|j>  =  f  s  j.(e)d£  =  /  S  +(e)dc  =  1  . 


(62') 


Now  m  =  iJj'iJj,  so  that  <m>  =  1  as  contended  above.  For  the  spectrum  of  m 
(not  an  analytic  signal)  we  find  by  convolution,  see  van  Vliet  [25], 


s;(f) 

m 


2  /  S  +(v)  S  ,  (v+  f)dv 

or  r 


=  2  /  S  , (v)  S  ,  (v  + f)dv 

-CO  *  +  *  + 


26(f)  +  4ctA(f/f0)aA  f_1u(f-f0) 

1  +  2aA  /A/h( f/f  )aA  f_1df 

f0 


where  a  small  ("noise  of  noise")  term  proportional  to  (aA)2  has  been 
neglected  both  in  the  numerator  and  in  the  denominator.  Again,  the  6(f) 
term  represents  the  DC  power  of  the  spectrum,  for  by  the  Wiener-Khintchine 
theorem, 
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S  *(f)  =2  f  e  <&m(t)Am(0)>  df 
Am 

—CO 

CO  oo 

=  2  /  e  iu)t  <m(t)m(0)>  dt  -  2  <n>  f  e  dt 

—CO  — .00 

=  S-(f)  -  2  6(f)  ,  (65) 

m 

where  we  noted  <m>2  =  1.  With  aA<<l,  (64)  and  (65)  give 

n  A  —  1 

SA'(f)  =  4  aA(f/fQ)aA  f\f  fQ  .  (66) 

Substituting  into  (61)  we  find  for  the  shot  noise  plus  1/f  noise  in  the 
flux  m(t) : 

SAm(w)  =  2  <m>  +  4  <m>2  aA(f /fQ)aAf _1  =  2  <m>  +  4  <m>2  aA/f.  (67) 

This  expression  ignores,  however,  the  incoherence  of  the  1/f  noise  con¬ 
tributions  from  various  carriers  in  the  flux.  If  n  is  the  concentration 
of  electrons  in  the  cathode  of  surface  area  F  and  X^  is  the  mean  free  path 
of  the  electrons,  the  number  of  independent  noise  contributions  is 
N  =  nFA^.  Hence  (67)  must  be  corrected  to  read 

S.  (w)  =  2  <m>  +  4  <m>2  aA/fN  .  (68) 

Am 

If  I  =  em,  then  for  the  current  noise, 

sAl(w)  =  2  e  <I>  +  4  <I>2  aA/fN  .  (69) 

Thus  we  have  a  Hooge  constant: 

ay  =  2aA  (70) 

For  a  computation  of  aA  for  various  infra-quanta,  see  Handel  [16]. 
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D.  Publications  and  Theses  Resulting  from  the  Contract 

1)  M.Sc.  Theses 


a)  R.R.J.  Schmidt:  "Various  Systems  for  the  Measurement  of  Noise 
from  20  Millihertz  to  30  Megahertz,"  University  of  Florida, 
March,  1981. 

b)  J.  Kilmer:  In  preparation. 

2)  Ph.D.  Dissertation 


H.  Mehta,  "Transport  Noise  Arising  from  Diffusion  and  Bulk  or 

Surface  Generation-Recombination,"  University  of  Florida,  January, 

1981. 

3)  Publications 

a)  "Theory  of  Transport  Noise  in  Semiconductors,  A  Survey," 

K.M.  van  Vliet  and  H.  Mehta,  Physica  Status  Solidi  (b) ,  to 
appear  July,  1981. 

b)  "Superstatistical  Emission  Noise,"  K.M.  van  Vliet,  P.H.  Handel 
and  A.  van  der  Ziel,  Physica  A,  in  press. 

c)  "Absence  of  Temperature  Fluctuations  in  1/f  Noise  Correlation 
Experiments,"  J.  Kilmer,  E.R.  Chenette,  K.M.  van  Vliet  and 
P.H.  Handel,  submitted  to  Phys.  Rev.  Letters. 

d)  Three  papers  are  in  preparation  on  the  work  done  by  H.  Mehta 
and  K.M.  van  Vliet. 


-  34  - 


E.  References 


t  1]  S.  Machlup,  "Scale  Invariance  Implies  1/f  Spectrum,"  Second  Inter¬ 
national  Conf.  1/f  Noise,  Orlando,  Florida,  1980,  Report  (Univer¬ 
sity  of  Florida),  p.  556. 

[  2]  P.H.  Handel  et  al.,  "Towards  a  More  General  Understanding  of  1/f 
Noise,"  submitted  to  Physics  Letters  . 

[  3]  R.F.  Vos  and  J.  Clarke,  Phys .  Rev.  B13,  556  (1976). 

f  4]  J.W.  Eberhard  and  P.M.  Horn,  Phys.  Rev.  Letters  39_,  643  (1977). 

[  5]  P.  Dutta,  J.W.  Eberhard,  and  P.H.  Horn,  Solid  State  Comm.  2]_,  1389 
(1978). 

[  6]  P.  Dutta,  P.  Dimon,  and  P.H.  Horn,  Phys.  Rev.  Lett.  4J3,  646  (1979). 

[  7]  K.M.  van  Vliet,  A.  van  der  Ziel,  and  R.R.  Schmidt,  J.  Appl.  Phys.  51, 
2947  (1980).  ““ 


t  8]  K.M.  van  Vliet  and  H.  Mehta, "Theory  of  Transport  Noise  in  Semiconductors 
A  Survey,"  Phys.  Stat.  Solidi  (b) ,  in  press. 

[9]  A.  van  der  Ziel,  "Flicker  Noise  in  Electronic  Devices,"  in  Advances 
in  Electronics  and  Electron  Physics  (L.  Martin  Ed.),  1979. 

[10]  H.  Mehta  and  K.M.  van  Vliet,  to  be  published. 


[11] 

[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 
[19] 


F.H.  Hooge  and  J.L.M.  Gaal,  Philips  Res.  Rep.  26_,  345  (1971). 
F.H.  Hooge  and  L.K.J.  van  Damme,  Phys.  Lett.  66A,  315  (1978). 
T.G.M.  Kleinpenning,  Physica  77,  78  (1974). 


R.P.  Jindal  and  A.  van  der  Ziel,  J.  Appl,  Physics,  in  press. 

P.H.  Handel,  Phys.  Rev.  Lett.  34,  1492,  1495  (1975). 

P.H.  Handel,  Phys.  Rev.  22A,  745  (1980):  also  Proc.  Second  Inter. 
Symposium  1/f  Noise,  Orlando,  Florida,  1980,  pp.  42,  56,  96,  525, 
and  550;  K.L.  Ngai,  ibid.,  p.  445. 


J.H.  Scofield,  D.H.  Darling  and  W.W.  Webb,  Phys.  Rev.  Letters, 


submitted 


M.B.  Weissman,  Phys.  Rev.  Letters,  submitted. 

See  [8]  and  also  Proc.  Second  Intern.  Conf.  1/f  Noise,  Orlando, 
Florida,  1980,  Report  (University  of  Florida),  p.  305. 


[20]  K.M.  van  Vliet  and  J.R.  Fassett,  in  Fluctuation  Phenomena  in  Solids 
(R.E.  Burgess,  Ed.),  pp.  267-359. 


[21]  G.G.  MacFarlane,  Proc.  Phys.  Soc.  (London)  B63,  807  (1950). 

[22]  R.E.  Burgess,  Proc.  Phys.  Soc.  (London)  B66,  334  (1953). 

[23]  K.M.  van  Vliet  and  E.R.  Chenette,  Physica  J33,  985  (1966). 

[24]  M.Lax  and  P.  Mengert,  Phys.  Chem.  Solids  14_,  248  (1960). 

[25]  K.M.  van  Vliet,  Physica  8J3  B  &  C,  52  (1976). 


30  dBV  FULLSCALE  -  10  dB/DIV  ZERO  COMP.  FILTER.  1 

AVERAGES.  64  32  THIRD-OCTAVE  SEGMENTS 


•30  dBV  FULLSCALE  -  10  dB/OlV  ZERO  COMP.  FILTER.  1 

AVERAGES.  64  32  THIRD-OCTAVE  SEGMENTS 


s 


S 

in 

i 


(9 

10 

l 


0  ‘00002 


0  ‘00001 


0  ‘0005 


0  ‘0052 


N 

0  ‘052 1  I 


> 

0  *595  Z 
Ul 

a 
a 
ui 
cc 

0  ‘2B2 


0  *521 


0  *E9 


3  'IE 


S3  ‘91 
S 


BO  NI  30VJ.10A 


ro 

ci 

•H 

Ui 


One-third  octave  display  of  collector  channel  noise  for  transistor  2N3945  111. 


Block  diagram  for  real  time  correlation  measurement  system 
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Two-stage  FET  low-noise  amplifier 
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Active  load  low-noise  preamplifier. 


VOLTAGE  IN  DB 


ca 

LO 

ca 

ca 

ca 

(9 

ca 

ca 

ca 

ca 

• 

• 

■ 

• 

■ 

• 

a 

a 

a 

• 

CO 

CO 

in 

CM 

in 

ca 

ca 

ca 

ca 

m 

CO 

CM 

oca 

CO 

in 

<a 

ca 

ca 

•—1 

CM 

in 

CM 

in 

ca 

ca 

*— • 

CM 

in 

ca 

FREQUENCY  IN  HZ 


Fig.  9 

Noise  spectrum  of  aluminum  film  at  1.071  mA 
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Fig.  12(b) 

Circuit  to  obtain  noise  spectra  and  coherence. 
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TRANSFER  FUNCTION 


Fig.  13 


Thermal  transfer  function. 
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wave  analyzer 


Three-transistor  heat-trans fer  measurement 


Thermal  transfer  function  for  the  two  "sensors 
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FIG.  19 

CORRELATION  (COHERENCE  FACTOR)  OF  THE  TWO  "SENSOR"  TRANSISTORS. 


0  dBV  FULLSCALE - 10  dB/DIV  ZERO  COMP.  FILTER.  25 

VERAGES.  4  32  THIRD-OCTAVE  SEGMENTS 
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Low-frequency  one-third  octave  display  of  base  channel  noise 
for  2N3945  //I  (see  also  Fig.  2). 


30  dBV  FULLSCALE  -  10  dB/DIV  ZERO  COMP.  FILTER*  25 

AVERAGES.  4  32  THIRD-OCTAVE  SEGMENTS 


2N3945  it  1  (see  also  Fig.  3)  . 
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Fig.  23 

One-third  octave  display  of  base  channel  noise  for  transistor  2N3945  // 1. 


Noise  spectra  for  embedded  sphere  Cn=0,  tv?0) 
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Fig.  29 

Noise  spectra  for  square  disc  and  square  bar 
(Physical  diffusion  noise  source) 
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